fid = fopen('diff_equ_1D_gas_phase_rec.dat');
data = textscan(fid,'%f %f %f %f %f','HeaderLines',15);
fclose(fid);

t = data{1};
Xs = data{2};
plot(t,Xs,'LineWidth',2);
set(gca,'FontSize',20);
xlabel('t(s)');
ylabel('Xs');
print('Xs','-dpng');
print('Xs','-depsc');